Transport Characterization of the Magnetic Anisotropy of (Ga,Mn)As 
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The rich magnetic anisotropy of compressively strained (Ga,Mn)As has attracted great interest 
recently. Here we discuss a sensitive method to visualize and quantify the individual components 
of the magnetic anisotropy using transport. A set of high resolution transport measurements is 
compiled into color coded resistance polar plots, which constitute a fingerprint of the symmetry 
components of the anisotropy. As a demonstration of the sensitivity of the method, we show that 
these typically reveal the presence of both the [110] and the [010] uniaxial magnetic anisotropy com- 
ponent in (Ga,Mn)As layers, even when most other techniques reveal only one of these components. 



The ferromagnetic semiconductor (Ga,Mn)As exhibits 
rich magnetic anisotropy behaviour. Experimental stud- 
ies of compressively strained (Ga,Mn)As based on direct 
magnetization measurements by SQUID magnetometry 
[H show a principally biaxial anisotropy in the sample 
plane for highly doped samples at 4 K. Additionally, 
uniaxial anisotropy components along [110] Q or [010] 
[1] have been reported. The [110] uniaxial anisotropy is 
widely accepted to be present in most (Ga,Mn)As layers 
as a secondary component at 4 K which gets stronger 
upon annealing, or increasing sample temperature. The 
primary biaxial anisotropy is expected from theory 0, Q 
and originates from the hole-mediated ferromagnetism 
and the strong spin orbit coupling, which links the mag- 
netic interactions to the lattice structure. This theory 
also provides a partial understanding of the emergence 
of a different easy axis for higher carrier concentration or 
temperature, in that it predicts a transition from biaxial 
along (100) to biaxial along (110). It does not however 
account for the presence of either uniaxial term. Never- 
theless, the existence of these uniaxial easy axes is ex- 
perimentally well established, and it remains of great 
interest, both as a practical matter in designing and 
controlling magnetization in experimental samples and 
as a fundamental tool for increasing our understanding 
of the magnetization behaviour in strong spin-orbit sys- 
tems, to develop tools to precisely investigate the mag- 
netic anisotropy in (Ga,Mn)As. 

Traditionally the anisotropy in these materials has 
been investigated by direct measurement of the projec- 
tion of the magnetization onto characteristic directions 
using SQUID or VSM. The advent of vector field mag- 
nets has recently opened up new possibilities for acquir- 
ing a detailed mapping of the anisotropy. We present one 
such method here, based on summarizing the results of 
transport measurements into color coded resistance polar 
plots (RPP) which act as fingerprints for the anisotropy 
of a given structure. Not only is this method faster than 
the traditional alternatives, but it is also more sensitive 
to certain secondary components of the anisotropy, in 
particular those with easy axes collinear to the primary 
anisotropy [§]. The technique thus often reveals the exis- 



tence of components which would be missed using stan- 
dard techniques. Moreover, the technique can be applied 
to study the anisotropy of layers by using macroscopic 
transport structures, or applied directly to device ele- 
ments. It can thus reveal any impacts of processing or 
the influence of small strain fields due to, for example, 
contacting. 

Here, we apply the RPP technique to verify that 
both the [010] and the [110] uniaxial anisotropy terms 
can be, and typically are, simultaneously present in 
compressively strained (Ga,Mn)As layers at 4 K. Color 
coded polar plots can be compiled from any measure- 
ment data that shows a response to the magnetization 
direction. Anisotropic magnetoresistance(AMR)p7), tun- 
neling anisotropic magnetoresistance (TAMR)[3] or in- 
plane Hall(IPH) measurements Q are typical examples in 
(Ga,Mn)As. For example, sending a current J through 
a (Ga,Mn)As Hall bar device along i9 = 0°, yields a 
sin(2'i^)-dependence of the transverse resistance Rxy on 
the magnetization direction i9 due to the IPH effect 0: 

We demonstrate the usefulness of the technique by ap- 
plying it to the characterization of a typical 20 nm thick 
(Ga,Mn)As layer grown on a GaAs (001) substrate by 
low-temperature molecular beam epitaxy. It is patterned 
into a 60 /im wide Hall bar oriented along the (110) crys- 
tal direction by optical lithography and chlorine assisted 
dry etching. Contacts are established through metal 
evaporation and lift off. 

Transport measurements are carried out in a magne- 
tocryostat fitted with a vector field magnet comprised of 
three mutually orthogonal pairs of Helmholtz coils that 
allow the application of a magnetic field of up to 300 mT 
in any direction. An angular set of IPH curves is ac- 
quired while sweeping the magnetic field along multiple 
directions in the sample plane. For each individual angle 
the magnetization state of the sample is first prepared 
by a strong negative magnetic field along cp. The field is 
then slowly brought down to zero while assuring that the 
field vector never deviates from the —(p direction. The 
IPH curve as a function of positive magnetic field in the 
cp direction is then acquired from zero to higher fields. 
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FIG. 1: (a) In-plane Hall measurement along Lp — 80° 
with marked first and second switching field, color scale and 
the corresponding section of a color coded resistance polar 
plot (inset), (b) Resistance polar plot from a full set of In- 
plane Hall measurements along every 3°. The 80°-section 
corresponding to (a) is marked by a dashed line. 

and these results are displayed in a RPP as in Fig. [T]3. 

Fig.[T^ shows a typical IPH measurement. After mag- 
netizing the sample along 80° at -300 mT, the magneti- 
zation relaxes towards the [010] easy axis, as the field is 
brought back to zero. This corresponds to the negative 
resistance state associated with an angle of = —135° 
between J and M. The first abrupt resistance change at 
the field i^ci corresponds to a 90° reorientation of M 
towards the other (Ga,Mn)As easy axis ~ [001]. A sec- 
ond reorientation of M towards [010] at i^c2 completes 
the magnetisation reversal. A set of such magnetic field 
scans along many angles, here every 3°, is compiled into 
a RPP with the magnetic field R along the radius and 
each scan at its angle (/:?. 

The inset of Fig. [T^ shows the 80°-segment of the full 
RPP of Fig.[T)3. The intensity encodes the normalized re- 
sistance value, where low and high denote the minimum 
and maximum resistance of the entire curve set, respec- 
tively. The positions of the switching fields in the polar 



plot and the symmetry of the pattern they form contain 
information on the underlying magnetic anisotropy. 

We now explain the observed pattern and the con- 
clusions that we may infer from it. As already noted, 
(Ga,Mn)As layers grown on a GaAs (001) substrate are 
compressively strained and exhibit a biaxial magnetic 
anisotropy (K\) with the two easy axes along [100] and 
[010] at 4 K. Additional smaller uniaxial anisotropy com- 
ponents collinear to one of the biaxial easy axes, for ex- 
ample [010] (i^ui) |3| or bisecting the biaxial easy axes, 
along [110] (i^n2)[2| have been previously reported. The 
energy £^ of a single (Ga,Mn)As domain with magnetiza- 
tion M pointing in 'd- and external magnetic field B. in 
(^-direction is thus phenomenologically expressed as: 

E = ^sin2(2^)+i^^i sin2(^)+i^^2sin2(^-45°)-Mi^cos(^-(/:)) 

^ (1) 
where the last term is the Zeeman energy. When an ex- 
ternal magnetic field is applied, the magnetization direc- 
tion follows a local energy minimum (coherent or Stoner- 
Wohlfarth rotation). However, if the energy gained by 
a magnetization reorientation to any other energy min- 
imum is larger than the respective domain wall (DW) 
nucleation/propagation energy £, a DW is nucleated and 
propagates through the structure, resulting in an abrupt 
magnetization reversal. Due to the mainly biaxial nature 
of the magnetic anisotropy the magnetization typically 
reverses in a double-step switching process through the 
nucleation and propagation of two 90° DW, as in Fig.[T^. 

We first discuss the regime, where DW nucle- 
ation/propagation dominates the magnetization reversal 
process, i.e. where e is much smaller than the crys- 
talline anisotropy. We start with the pure biaxial mag- 
netic anisotropy, i.e., with K^i = i^n2 = 0. During the 
double-step switching process, the magnetostatic energy 
minima remain to a good approximation along the bi- 
axial easy axes, whose difference in energy is given by 
the Zeeman term in Eq. [H When the energy gained 
through a 90° magnetization reorientation is larger than 
£900, the nucleation and propagation energy of a 90°- 
DW, a thermally activated switching event becomes pos- 
sible, which on the timescale of our measurement, re- 
sults in an immediate switching event. The polar plot in 
Fig. [2^ shows the characteristic square pattern formed 
by these switching fields 0, Fig. [2]3 the correspond- 
ing calculated RPP. The switching field positions (thick 
black lines) can be expressed in cartesian coordinates 
{x = Rc cos i^\y = iJc sin ip) as given in Tab. [D The 
diagonals (grey) of the i^ci -square represent the biaxial 
easy axes of the material(here along 0° and 90°). The 
diagonal's length is equal to . The dashed lines rep- 
resent the hard magnetic axes. The arrows illustrate the 
direction of the magnetization and their color the corre- 
sponding resistance state of the respective section. 
An additional small uniaxial anisotropy K^i along one 
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FIG. 2: (a) Switching field positions (thick sohd hnes) in a 
polar plot for a biaxial material with easy axes along [100] and 
[010] (gray). The magnetization direction in each region of the 
plot is indicated by arrows (red/black: high/low resistance) 
and the hard axes by dashed lines, (b-d) Simulated resistance 
polar plots for a biaxial material with easy axes along the [100] 
and [010] crystal directions (b) and the same material with 
an additional uniaxial anisotropy along [010] (c) or [110](d). 
Color scale of the resistance as in Fig. [1] e denotes the 90°- 
DW nucleation/propagation energy. 



Anisotropy components 


ifc-positions in polar plot 


(100) 
(100) + [010] 

(100) + [TlO] 


y — M ^ M 

_ f+x± V2^cos(45° - V2) 
^ ~ 1 -X db cos(45° + 6/2) 



TABLE I: Positions of the switching fields in a polar plot 
given in cartesian coordinates (x = He cos (p;y = He sin ip). 
6/2 is the angle between global and biaxial easy axis as in 
Fig. 01 (see text). 



of the biaxial easy axes (here along 90°) alters this pat- 
tern as shown in Fig. The symmetry between the 
two biaxial easy axes is lifted, since one of them is par- 
allel {Belle) and one perpendicular (Bellh) to the easy 
axis of the uniaxial component. The switching field po- 
sitions in Tab. [Jean be derived as discussed above [6|. 
90°-switches away from (towards) the Belle axis occur 
now at higher (lower) magnetic fields as compared to 
the pure biaxial anisotropy. Typical steps in the Hd- 
pattern emerge (Fig. [2]3, e.g. along 45°). The strength 
of the uniaxial anisotropy can be determined from the 
separation between Hci and Hc2 along the BeUh 

axis. Another characteristic feature is the "open corner" 



of the Hci- pattern along the Belle axis. A 180° magne- 
tization reversal becomes energetically more favorable in 
this angular region, than two successive 90° events [gI, [H . 
Since the isotropic magnetoresistance of typical samples 
is relatively small, two magnetization directions differing 
by 180° are not distinguishable on the scale considered 
here, and have the same color in the RPP. 

If on the other hand, apart from the main biaxial 
anisotropy (i^i,(100)), a uniaxial anisotropy Ku2 with 
easy axis along [110] (here 135°) is also present and 
assuming that the DW nucleation/propagation energy 
depends on the angle Ai9 between the two domains as 
£Ai9 = ^90° (1 — cos(A'i^)), the calculation gives the RPP 
of Fig. [2]i. The uniaxial anisotropy component flattens 
the energy surface and shifts the positions of the biaxial 
energy minima by | = ^arcsin(^^^) towards the uni- 
axial easy axisQ. The resulting "global easy axes" ap- 
pear as diagonals of a characteristic rectangle in the RPP, 
where the uniaxial easy axis is along the median line of 
the longer edge. The positions of the switching events 
(Tab. [Ij) can be derived equating e with the difference 
in Zeeman energy between the minima. The diagonal's 
length is again equal to . 

We demonstrate the usefulness of the technique by ap- 
plying it to the characterization of the typical (Ga,Mn) As 
sample discusses in Fig. [TJ The RPP in Fig. [TJd shows an 
i^ci-pattern with both an elongation in [110]-direction 
and steps along the hard axes. Both uniaxial anisotropy 
components are thus clearly present. From the length of 
the diagonals we estimate ~ 8 mT. The step height 
gives ~1 mT. From the rectangle side-ratio we ob- 
tain (5-8°. Ku2 is thus -15% of i^i. ^ -100 mT 
can be estimated from the asymptotic behavior of the 
magnetization towards the hard axes at higher fields and 
M — 37,000 A/m is known from SQUID measurements. 

The anisotropy components and e can differ from sam- 
ple to sample. Neglecting coherent rotation is typically 
a good model for the first switching fields whereas 
Hc2 is strongly influenced by magnetization rotation es- 
pecially along the hard axes, as can be seen in Fig. [TJd. 
Sharp switching events are easily identified by high con- 
trast, whereas gradual magnetization rotations show up 
as smooth color transitions. 

In conclusion we have shown that RPP compiled from 
high resolution IPH measurements constitute a finger- 
print of the (Ga,Mn)As anisotropy. They allow both 
qualitative and quantitative statements about the sym- 
metry components of the magnetic anisotropy and the 
DW nucleation/propagation energy. The same technique 
is equally applicable to any transport phenomena which 
produces a response to the orientation of the magnetiza- 
tion, such as AMR, or TAMR. 

We applied the method to several typical MBE lay- 
ers from our lab [13] to confirm that biaxial and both 
the [110] and the [010] uniaxial anisotropy term are 
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present in all typical compressively strained (Ga,Mn)As 
layers at 4 K with a relative strength of the order of 
Ki : Ku2 • Kui - 100 : 10 : 1. Indeed ah (Ga,Mn)As 
layers investigated show both these uniaxial components, 
including layers where the [010] component could not be 
identified in SQUID measurements. Moreover the appli- 
cation of our fingerprint method to previously published 
data in the literature shows that in all cases where suf- 
ficient data is available, both uniaxial components are 
present. 
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